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We study near a Feshbach resonance, as a function of the mass ratio, the fermion-dimer scattering 
amplitude in fermionic mixtures of two fermion species. When masses are equal the physical situation 
is known to be quite simple. We show that, when the mass ratio is increased, the situation becomes 
much more complex. For the s-wave contribution we obtain an analytical solution in the asymptotic 
limit of very large mass ratio. In this regime the s-wave scattering amplitude displays a large number 
of zeros, essentially linked to the known large value of the fermion-dimer scattering length in this 
regime. We find by an exact numerical calculation that a zero is still present for a mass ratio of 15. 
For the p-wave contribution we make our study below the mass ratio of 8.17, where a fermion-dimer 
bound state appears. We find that a strong p-wave resonance is present at low energy, due to a 
virtual bound state, in the fermion-dimer system, which is a forerunner of the real bound state. 
This resonance becomes prominent in the mass ratio range around the one corresponding to the 
40 K - 6 Li mixtures, much studied experimentally. This resonance should affect a number of physical 
properties. These include the equation of state of unbalanced mixtures at very low temperature but 
also the equation of state of balanced mixtures at moderate or high temperature. The frequency and 
the damping of collective modes should also provide a convenient way to evidence this resonance. 
Finally it should be possible to modify the effective mass of one of the fermionic species by making 
use of an optical lattice. This would allow to study the strong dependence of the resonance as a 
function of the mass ratio of the two fermionic elements. In particular one could check if the virtual 
bound state is relevant for the instabilities of these mixtures. 

PACS numbers: 03.75.Ss, 03.65.Nk, 34.50.Cx, 67.85. Lm 

I. INTRODUCTION 

One of the most fascinating aspects of the physics of ultracold gases pQ is the ability, in appropriate situations, to 
describe the interaction between non identical fermionic atoms in terms of a single parameter, namely the scattering 
length a for two different atoms. For example in fermionic gases involving only atoms in the two lowest hyperfine 
states of a single element, such as 6 Li or 40 K, the interaction is fully described by the scattering length between atoms 
belonging to these two different hyperfine states [1] . Moreover in the presence of a Feshbach resonance this scattering 
length can be experimentally modified almost at will by merely changing the applied magnetic field. The case of wide 
Feshbach resonances is particularly convenient since all the physics related to the origin of the resonance, namely 
the existence of a closed channel, is irrelevant for most purposes. This has allowed, in the case of two fermionic 
atomic species, the experimental realization of the BEC-BCS crossover where the scattering length goes from small 
negative values to small positive values through the resonance, where it diverges. For small negative values of the 
scattering length two different atoms have a weak effective attraction, which for high enough density and small 
enough temperature, leads to the formation of Cooper pairs analogous to the ones appearing in superconductors and 
corresponding superfluid properties. On the other hand for small positive values of the scattering length a > 0, two 
fermionic atoms of different species form a bound state resulting in the appearance of molecules, or dimers, with 
binding energy E^. These dimers behave statistically as bosons, and again for high enough density and small enough 
temperature a Bose-Einstein condensate forms with superfluid properties. The BEC-BCS crossover occurs when one 
goes continuously from the Bose condensate to the BCS superfluid by merely varying the scattering length through 
the change of magnetic field. 

The simplicity of the description of the interaction extends to the BEC limit, when one considers dimer-dimer 
interaction, while this is a priori a more complex situation since dimers have an internal structure. Nevertheless the 
BEC limit is quite simple since the composite bosons one deals with are dilute. In this case again the interaction 
between these bosons is fully characterized by their scattering length 04. This is because in the low temperature 
regime corresponding to the existence of the Bose-Einstein condensate, the kinetic energies of the dimers will be quite 
small compared to their binding energy, and the structure of these composite bosons will be irrelevant. This scattering 
length 04 has been obtained 0[3], in terms of the scattering length a for fcrmions, as 04 ~ 0.6 a when the different 
fermions have equal masses. 

On the other hand when one goes away from the BEC limit toward the unitarity limit, it is no longer true that the 
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composite nature of the bosons is irrelevant. One might accordingly expect a more complex physics for this molecular 
gas. Remarkably the T = equation of state on this BEC side is nevertheless very well described [1] in a wide domain 
of density by the Lee-Huang- Yang [J] equation of state for these composite bosons [5] , making again merely use of the 
dimer-dimer scattering length 04 between these bosons. This is as if these bosons were elementary. This is somewhat 
surprising. 

More precisely, assuming for example two fermions species with equal densities n — fc^/67r 2 , a typical fermion 
wavevector is of order kp, and for kpa ~ 1 the typical kinetic energy of a fermion is now of the same order as the 
dimer binding energy which is Ei, = l/(2/^a 2 ) with \i being the reduced mass. In this case one has to consider the 
whole energy dependent scattering amplitude and not only its zero energy limit given by the scattering length. This 
is valid for dimer-dimer scattering as well as for single fermion-dimer scattering. In this last case these single fermions 
may arise because of thermal breaking of dimers, or because we consider unbalanced fcrmionic mixtures. Actually for 
kpa ~ 1 these single fermions arise also because we may no longer consider isolated dimers since they are too close, 
and we have actually to deal with a complicated many-body problem. However our point is that even the simplest 
minded description has to consider the energy dependence of scattering processes. 

Actually it is possible to understand qualitatively the wide range of validity of the Lee-Huang- Yang equation of 
state in the following way. One assumes first that only the s-wave scattering properties are relevant, since higher 
partial wave give a quite small contribution at low energy. Then one supposes that the s-wave phase shift So(k) stays 
small enough so that its low energy expansion So{k) ~ — ka remains valid, where a is the scattering length relevant 
for the considered process, dimer-dimer or fermion-dimer scattering. In this case the standard low energy expression 
for the scattering amplitude f(k) = — 1/(5 _1 + ik) remains valid. This is what happens when the different fermions 
have equal masses. In this case the scattering amplitude has a fairly weak variation when k is going from zero to the 
maximum value corresponding to dimer breaking, so the situation is not much different from evaluating the scattering 
amplitude at k = 0. 

The purpose of the present paper is to show that this simplicity is linked to the fact that the two atomic species 
have equal masses = to^ (we follow the standard use of referring to the two different atomic species as f and i 
atoms, even if there is no relation between this notation and the physical spin of these atoms). It disappears when 
the masses are quite different. In this case, even though the elementary scattering process between different fermions 
is very simple to describe, nevertheless the mere complexity due to the existence of a dimer leads to complexity in the 
involved scattering processes. In the present paper we will consider only the simplest situation which displays such a 
feature. We will deal with the fermion-dimer scattering, not the dimer-dimer scattering. The lonely fermion has an f 
spin and scatters with a f-i dimer. We will furthermore assume that the available kinetic energy is small enough so 
the dimer is not broken by the collision. In other words we will deal only with elastic scattering. 

In a first part of the paper we will consider s-wave scattering, since it is the dominant one at very low energy and is 
anyway expected to be an essential component of the scattering process. This is the natural extension of the case where 
the scattering is fully described by the scattering length. Actually, although we have performed numerical calculations 
for any value of the mass ratio to^/to,^, we deal extensively with the limit where this mass ratio is large. Indeed in 
this case we will be able to obtain an analytical solution of the problem, which has the advantage of displaying quite 
explicitly the behaviour resulting from the mass difference between the two fermion species. Naturally since we have 
already explored in a preceding work [5] the behaviour of the scattering length in this limit, this paper can also be 
seen as a natural extension of our earlier work. But it has the major interest of displaying a strong structure in the 
s-wave scattering amplitude, which is actually essentially linked to the increase of the scattering length with the mass 
ratio [SHE]- 

Then in the second part of our paper we will turn to p-wave scattering. Although it does not contribute in the zero 
energy limit, this partial wave turns out to be quite important. This is physically linked to the appearance of new 
bound states which appear at the level of the three-body problem when the mass ratio increases. These bound states 
have been studied in great details by Kartavtsev and Malykh [Sj, for the present situation of two identical fermions 
with mass and a different fermion with mass mj,. They are not present in the s-wave channel and appear only 
for angular momenta £ > 1. Among them are the well-known Efimov |10j states, which are quite remarkable for their 
spectrum structure with binding energy going up to infinity when the range of the interaction goes to zero. For £ = 1 
the Efimov states appear only [HUH] for a mass ratio to f /to I— 13.6. However there are already bound states for 
smaller mass ratios, which are not Efimov states. A first bound state [5] appears for mass ratio to t /m 1= 8.17. 

Actually very important contributions in the fermion-dimer scattering properties arise even for quite lower mass 
ratios. The basic reason is that, even before the mass ratio threshold for bound states is reached, virtual bound states 
are already present. For nonzero angular momenta, in particular £ = 1, the centrifugal barrier acts to inhibit their 
decay and provide them with fairly long lifetime. Hence they are not much different from real bound states. Moreover, 
since their energy is positive, they give rise to resonances, which are strong when they are located at low energy. In 
particular the p-wave contribution becomes very strong, due to quasi-resonance, when the mass ratio approaches 8.17. 
However we find that these effects are already important for the mass ratio m j" /to 1= 6.64 corresponding to the 



3 



mixtures of 6 Li and 40 K, which is probably presently the most investigated experimentallyJT^HH] and which is of very 
high current interest [15 - 23 . Actually this point has already been noted, for the specific case of the 6 Li- 40 K mixtures, 
by Levinsen et al [53]. This was mostly done for the case of narrow Feshbach resonances, which are by far the most 
frequently found in these mixtures. They were also mostly interested in lower dimension physics. 

Here we will consider only the case of wide Feshbach resonance in 3D, because of its theoretical interest due to 
its simplicity. In the following we will assume that we are near such a Feshbach resonance and ignore any stability 
problem. Fortunately such a Feshbach resonance, which is fairly broad and reasonably stable, has been identified 
experimentally near 155G [25] in 6 Li- K mixtures. We will consider in detail the mass ratio dependence of the p- 
wave contribution. This is not only of theoretical interest since we will show that, by making use of optical lattices, it is 
possible to vary experimentally the effective mass ratio of 6 Li- 40 K mixtures. This would allow to study experimentally 
the effect of the p-wave contribution when it is quite strong, or conversely to eliminate it to a large extent which could 
be in particular be quite useful if this contribution turned out to be experimentally detrimental. We note indeed 
that these virtual bound states could play a very important role in the three-body decay and in the stability of these 
mixtures, since they correspond to physical situations where three atoms stay close together for a fairly long time, 
increasing in this way the probability of three-body decay processes. In this way optical lattices could also be used 
to study and possibly remove instability problems. 

In this p-wave study we will consider mass ratios below the value 8.17, corresponding to the appearance of the 
first p-wave bound state, since clearly beyond this threshold the physical situation will get more complicated. This 
justifies also that we consider only the s-wave and p-wave components of the full scattering amplitude, since the higher 
angular momentum components are expected to give very small contributions to the total scattering in the regime we 
investigate. 

In practice the next section [TT] is devoted to the study of the s-wave contribution to the scattering amplitude. But 
in the first subsection we will present the basic general equations which will be used also for the p-wave contribution. 
Then in the rest of this next section [IT] we will study analytically in details the case where the mass ratio m f /m J, 
is very large. We will conclude this section by considering numerically cases where this mass ratio corresponds to 
situations which can be reached experimentally. Then in the last section |Hl| bcfore our conclusion we study numerically 
the p-wave contribution. 



II. S - WAVE SCATTERING 



A. Basic equation 

In the same way as in our earlier work [5] we use the integral equation formulation of this dimer-fermion scattering 
problem first worked out by Skorniakov and Ter-Martirosianj2(J]. It is convenient to derive it rapidly in the case where 
the fermion masses are different (27) by writing directly [3] the integral equation satisfied by the full dimcr-fcrmion 
scattering vertex T 3 (pi,p 2 ; P). Here P is the (conserved) momentum-energy of the particles P = {P, E}, with P the 
total momentum and E the total energy of incoming particles. p\ is the momentum-energy of the incoming f fermion 
and P2 its outgoing value. Since we want the dimer-fermion scattering amplitude, we can take the center of mass 
referential for which P — 0. If we call ko the outgoing momentum of the f fermion, i.e. P2 = ko, the total energy E 
is the sum of the binding energy of the dimer and of the kinetic energies of the outgoing particles: 

1 kl 

2/ia 2 2jUt ^ ^ 

where fi — m^m^/(m^ + is the dimer reduced mass, /i-r is the atom-dimer reduced mass fix — m^M/(m^ + M) = 
m^(m^+mi)/(2m^+m^). Here M — m^+m^ is the mass of the dimer. The integral equation satisfied by T 3 (pi,p 2 ; P) 

is ®-. 

T 3 (p uP2 ;P) = -G^P - Pl -p 2 )-J2 G^P - Pl - q)G t (q) T 2 (P - q) T 3 (q,p 2 - P), (2) 

where ^ = i J dq<if2/(27r) 4 , with q = {q, £1}. Here T 2 (P) is the dimer propagator: 

y rp\ _ \ rg\ 

M a- 1 - y/2^{V 2 /2M - E) 

while G^^(p) are the single fermion j" or I propagators. Just as in the equal mass case [3] one can show that, in the 
right-hand side of Eq.Q, only the on-the-shell value of T 3 (q,p 2 ; P) with respect to the variable q is needed. Hence 
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we may restrict ourselves to consider only this on-the-shell value everywhere in Eq.Q. Finally we have also to take 
the on-the-shell value for p2 = {p 2 , W2} to obtain the scattering amplitude, which means ui2 = k\j2m^. We are thus 
led to consider more specifically T 3 ({k, k 2 /2m^}, {k , k 2 /2m^}; {0, E}) and the scattering amplitude will be obtained 
by taking |k| = ko. Hence we introduce, instead of T^(pi,p2', P), a function d3(k, ko) which is directly related to the 
scattering amplitude, as we will see below [2*S] : 



fl3(k ' ko) "2^ 




i + ^!( fc 2_ fco2) 



T 3 ({k, fc 2 /2m T }, {k , k 2 /2m t }; {0, E}) (4) 



In this part we are only interested in the s-wave component of the scattering amplitude. This component is obtained 
from the above T3 by averaging it over the angle between the incoming momentum k and the outgoing one ko, and we 
note T3 the resulting angular average. Noting merely as(k) the corresponding average on 03 (k, k ) (naturally 03 (fc) 
still depends on ko = |ko|), we see that it is simply given by Eq.Q where T3 has been replaced by T 3 . Obviously 
when we take the angular average of integral equation Eq.Q, we are left with an integral equation involving only T3, 
or equivalently 03. 

It is more convenient to write the resulting integral equation by making use of reduced units, obtained by setting 
k = k/a, k — ko/a, q = q/a and a^{k)/a = 5 3 (fc). We obtain: 



3 (fc) _ 1 1 - Rk 2 , + k + K + R ' k h 



1 + ^1 + R ( fc 2 _ j.2) 2R'kh l-Rkl + k? + kl- R>kk 1 ' 

1 f°° d _ fa 3 (q) 1 1 - Rkl + k 2 + f + R'kq 
ttR' J q q 2 -k 2 - id kq ™ 1 - Rk 2 + k 2 + q 2 - R'kq 

where we have used for the involved mass ratios the simpler notations R = n/ ht, R' = 2[i/m±. In writing this 
equation we have used the fact that we consider only the case where the incoming dimer and fermion have kinetic 
energies low enough so that the dimer can not be broken in the scattering process. This implies that E < 0, or i?fc 2 < 1 
in our reduced units. This makes all the arguments in the logarithms positive. Otherwise we should have included 
— iS contributions in these arguments. In the case where ko — > one checks easily that this equation reduces to the 
one studied in [6 , which gives the scattering length for the same process. It is clear that B 3 (fe) is a complex quantity. 
From Eq.([5]) one can write the coupled integral equations for its real and imaginary parts, which are convenient for 
the numerical solution of Eq. ^ . In the case where the masses are equal = one can check that the resulting 
equations are identical to the ones obtained by Skorniakov and Ter-Martirosian_2B] for this specific case. 

Finally it is convenient to consider the relation between 03 (k, ko) and the wave function V'k(r) for the dimer-fermion 
scattering problem, as it is found in [26] and [28]. For large dimer-fermion distance r it behaves as: 

^ k (r)^e 4k - r + / fc (0)^ (6) 

where 9 is the angle between the incoming fermion momentum k and its outgoing value ko (with |k| = |ko|), and 
fk{0) is the scattering amplitude. Its Fourier transform reads: 

^ k (q) = J dr e -^ r Vk(r) = (2ir) 3 6(k - q) - -fr^Z^ (7) 

Note that the second term has a sign opposite to the one found in [26] . just because our definition for 03 has an 
opposite sign. The scattering amplitude is obtained from fk (&) — ~ 03 (q, ko) (with |q| = |k |). 

When one takes the angular average of Eq.(|6]), the resulting s-wave component tp^.(r) of the wave function is given 
by: 

^° k (r) = ^ r sm(kr + S ) (8) 
where 5o(k) is the s-wave phase shift, linked to s-wave component fo(k) of the scattering amplitude by: 

2i8 _ 1 

AW - -^r- (9) 
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B. Very light mass 



Let us go now to the interesting regime where the mass of the lonely fermion is very light (or equivalently the 
mass of the two identical fermions is very heavy). In this case parameter R is near zero and R is near 2. In 
our preceding paper [B], where we had fco = 0, we have carefully expanded Eq.|5j near this limit. Then by a Fourier 
transform we have converted the resulting equation into a differential equation with respect to a variable which is 
called r below. We have studied this last equation, matching in particular the behaviour for the small and the large 
values of the variable r by appropriate consideration of the boundary conditions. However the important range, in 
order to obtain the scattering length, corresponds to the large f values and the proper behaviour in this range can also 
be obtained from physical considerations, as we will see below. We can also proceed by continuity, requiring that in 
the present case the large variable r behaviour reduces to the one we have found in the limit fco — > 0. More specifically 
we have seen in this preceding work that the careful expansion of the right-hand side of Eq.|5]) was useful only to 
satisfy the perfect matching with the boundary conditions for the small values of the variable f. On the other hand 
only the expansion of the left-hand side of Eq. ^ was necessary in order to obtain the large f information necessary 
to obtain the scattering length. Here we will take advantage of this fact to avoid the painfull task of expanding the 
right-hand side of Eq.([5]) and rather rely on the above arguments to obtain the required large variable f behaviour. 
Hence we will set R and R' to their limiting value, namely R — and R' — 2 in the right-hand side. In the left-hand 
side, following our preceding work, we will also set 1 + R (k 2 — fcp) ~ 1 because large values of k 2 — kg are irrelevant 
in our problem. However we will naturally keep the overall factor R in front of the left-hand side. This is the only 
difference with taking the full = limit. 

We rewrite Eq.([5]) in this limit by making use, instead of 03(9), of the more convenient function F(q) defined by: 



q 2 — k 2 — iS 



m - ... do) 



We note that, since from Eq.Q as(q) is an even function of q, F(q) is an odd function of q. Then Eq.([5]) becomes 

l_ 1 + (fc + fc ) 2 
4fc n i + (fc-fc ) 2 



^-IDFik) - ^mi±##S (ii) 



^ ^ dq F(q) In ) + g + g ^ = ±- f° dq F(q) In (l + (k - qf) 
2tt Jo 1 + (« - q) 2 2tt J^^ 



where in the last equality we have made use of the odd parity of F(q) and, as in [B], we have set: 

9 FL m . „. 

e 2 = V " — (12) 

We take now the Fourier transform of this equation, introducing: 



F(F) = / dqF(q) exp(-iqf) (13) 

J — 00 

and making use of: 

r°° e -f 

dq In (l + q 2 ) exp(-iqf) = -2tt^- (14) 

where we have restricted ourselves to the case f > 0. We obtain: 

^ ( £ + ^ w = e 4- ( m + i.^-) (15) 

Since the e — solution: 



df 2 v w f V fco 



- . sm(fc r) 

F a{r) = -m — f (16) 

fco 

satisfies (d 2 /df 2 + fco)F (f) = 0, it is convenient to set: 

F(f) = F Q (f) + iirg(f) (17) 
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to obtain the simpler equation: 



df 2 



9(r) 



-9(r) 



(18) 



The fact that we find a second order differential equation is directly related to our lowest order expansion of Eq. ^ 
in mi/m^. Going to higher order in this expansion would result in a higher order differential equation 

At this stage it is useful to notice the relation between g(f) and the scattering wavefunction Eq.([8]). Taking in 
Eq.([7| the angular average over the direction of k, performing the angular integration over r, setting ]k| = |k | and 
making use of the definition Eq.(lO), we obtain: 



° 2n 2 
4-7T / dr r sm{qr)\j/l (r) = — — [S(q 
/o *0 



k )-6(q + k )]-47rF(q) 



(19) 



where we have added in the right-hand side the term S(q + kg) which is usely zero since q — |q| > and fco = |ko| > 0. 
However we have extended above the range of variation of q to negative values and made use of the fact that the 
corresponding extension of F(q) is odd. The added term is necessary to satisfy in Eq.(19) this parity property of 
F(q). Going to reduced variables and taking the Fourier transform of Eq.(19) for f > we obtain with Eq.(16|: 



e iS ° 

F(f) = F (f) + itt-j— sin(fc r + <*„) 
fc 



Hence from Eq.(17) we have for large r the very simple relation: 



JSr, 



k 



sin(fc r -I- Sq) 



(20) 



(21) 



This provides us with the boundary condition necessary to solve the second order differential equation Eq.(18), which 
is accordingly an effective Schrodinger equa tion for our scattering problem. 



Except for the k 2 term this equation Eq.(18) is the same as the one we have found in our preceding work [5] and 
we solve it by following the same procedure. Performing the change of variable r = 21n(2/z), implying z — 2e~ r / 2 , 
we find that g(z) = g(f{z)) satisfies: 



d 2 g{z) , dg{z) 



dz 2 



+ z ■ 



dz 



4k 2 g(z) 



L{z) 



g(z) = 



(22) 



with L(z) = 2e 2 ln(2/z) 



We notice again that, in the range of variable which is of interest, that is large 



f, corresponding to small z, L{z) is a quite slowly varying function of z. This allows to treat it as a constant 
L{z) ~ L. As discussed in details in [5] this approximate treatment is increasingly accurate when the mass ratio 
m^/m-f, i.e. e 2 , becomes increasingly small, which is precisely the limit we are interested in. We have found [ 6] t hat 
the important range for matching the solution of Eq.([T8]) to the asymptotic behaviour, in the present case Eq.(21 1, is 
f ~ fo = ln(l/e 2 ). Hence we may take L — e 2 ln(l/e 2 ). 



Then the further change of variable z — y L x transforms Eq. ( 22 1 into: 



d g(x) dg(x) 



dx 2 



dx 



4fc g{x) - xg(x) = 



(23) 



for g(x) = g{z). The solution of this equation which matches, as we will see shortly, the asymptotic behaviour Eq.(21 ) 
is the well known Bessel functions K 2i j :o (x) (or equivalently K_ 2i k g (x)), which goes continuously to the solution 
K (x) that we have found in [5] for the case fco = 0. A convenient integral representation for this function is |29) : 



dtt 



-x cosh t 



cos(2fc i) 



From this representation one obtains the small x behaviour: 



2 ikn 



T(2ik ) 



(24) 



(25) 



where T(x) is the Euler function. From this result one recovers in particular, for fc — > 0, the known small x asymptotic 
behaviour Kq(x) ~ ln(2/a;) — C where C is the Euler constant. 
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Since 2/x = \fhe r l 2 we see that, for large r, our solution has indeed the required asymptotic behaviour Eq.(21|. 
This is naturally within an unimportant constant prefactor since the homogeneous differential equation Eq. ( 23 1 does 
not allow to find this prefactor. But, as soon as S is found, this prefactor is given explicitly by Eq.(21 1. Comparing 
Eq.(25| with Eq.(2f ) in order to identify the phase shift, we find easily: 



S (k ) = k InL + Arg \iT(2ik )] = k Q In [e 2 ln(l/e 2 )] + Arg \iT(2ik )] 



(26) 



In particular in the limit k — > we can make use of T(2ik ) ~ l/(2ik ) — C, implying Arg [ir(2ifc )] ~ ~2k C. 
Comparing to the limiting behaviour 5o(ka) ~ — a^kg we recover our result [6j: 



53S ^ = _l nL + 2C = ln r ^-lnln^ + 2C7 



(27) 



for the fermion-dimer scattering length in the considered limit where the mass ratio m^/m^ is large. 

Once the phase shift is known we obtain the scattering amplitude by merely making use of Eq. ([9| . This is obviously 
the simpler way. However let us briefly show that we recover this result by calculating directly 03 (q, kg) from our 
solution K 2ikg (x) of Eq.(23). From section II A we have: 



/o(fco) = = -lim f _, So a 3 (q) = -21im f _, So (q ~ k Q ) F(q) 



(28) 



the last equality resulting from Eq.pO|. In the Fourier transform of the decomposition Eq.( 17) we see that F (q) does 
not contribute to the right-hand side of Eq.(28) since, as we have seen in Eq.(19): 



Fo(q) = XTT [S(q-~ko)-S(q + k )} 
ZKq 



so that, in this right-hand side, we are left only with the contribution of g 

df 



foih 



-21im 



ko) 



2tt' 



(r) exp(igr) 



(29) 



(30) 



We have seen that the solution of Eq. ( 23 1 is proportional to K 2i g (x) . Going back to the r variable this translates 
into: 



7(r) = a{k )K 2i - kQ 



2e~ f / 2 



(31) 



where, as mentionned above, the prefactor a(ko) can not be obtained from the homogeneous equation Eq.(23) and 
has to be determined from the asymptotic behaviour Eq.(21) for g(f). From the asymptotic behaviour Eq.(25| we 
have: 



K. 



2e- f / 2 
IT 



|r(2t* )| sin(fc f + ^ ) 



where 5q is given by Eq.(26). Comparing Eq.(|31[) with Eq.(21) we find: 



a(k ) 



o iS _ 



1 



fc |r(2 l fc )| 



The Fourier transform of g(r) is then: 

rao df 



2tt 



(r) exp(igr) 



dr 
2^ 



g(r) exp(igr) - (q -> -q) 



(32) 



(33) 



(34) 



where we have made use of the odd parity of g(r) directly linked to the odd parity of F(q). Going back to the variable 
x = 2e~ 7 'I 2 1 ' \[~L we have: 



00 dr a(k ) 
~{r) expiiqr) — 



2tt- 



2tt 
a(ko) 



dr exp(iqr) K 2fko 

v -2iq 



2e- f ' 2 
y/L 



(35) 



dxx-^K^x) 




FIG. 1: £ = phase shift So(ko) as a function of fco = fco a, given by Eq.(26l (dashed line) and its asymptotic approximation 
given by Eq.(41l (full line) for a) m\/m\, = 10 and b) m-f/mi = 100. 
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where in the last step we have replaced the upper bound 2/\/Z by oo, since in our case L = e 2 ln(l/e 2 ) <C 1 and, 
as it is clear from Eq.(24|, K 2i k Q {x) decreases exponentially rapidly for large x. The last integral can be performed 
analytically [22], leading to: 



~ g(f) exp(igf) = ^-L~*T (i(k - q)) T (-i(k + q)) - (q -> -q) 



(36) 



Since T(x) ~ 1/x for x -> 0, we find from Eq.Q and Eq.(|36l 

a(fco) 



/o(fco) = 



L lfco r(2ifc ) + 2r lfen r(-2ifc ) 



However we have from Eq.(26|: 



sm(<5o(&o)) — £ p^+i 2|F(-2zfc )| 



Taking Eq.(33) into account this leads to: 



/o(fco) = e = 

fc 



(37) 



(38) 



(39) 



which, when we go back to unreduced units, is identical to Eq.([9]) as it should be. This shows that our calculations 
and approximations in the very light mass limit are fully consistent. 



C. Discussion 



Let us now consider the implications of our result Eq.(26) for the phase shift and the corresponding scattering 



amplitude. Although the Euler T function is easily obtained numerically, we have more insight if we make use of an 
excellent approximation. Since for real x we have: 



Arg [ir(ix)] = Arg [ixT(ix)] = Arg [r(l + ix)} = Im [lnT(l + ix)] 
we may use the excellent Stirling formula lnT(z) ~ (z — 1/2) hxz — z + (1/2) ln(27r) which leads to: 



(40) 



5 (fc ) = -k Q a 3 + - arctan(2fc ) + 2(C - l)fc + fc ln(l + 4fcg) 



(41) 



If we were to correct the small imperfection of the Stirling formula for low z by replacing 2C ~ 1.154 by 1, this would 
allow Eq.(41| to give the exact result livaj. _ ¥0 6o(ko) = —koa^ with S3 given by Eq.(27), but the result would not be 
accurate for higher values of k$. Our original result Eq.(26l and its simplified form Eq.(41) are displayed in Figfljfor 
m-f/m^ = 10 and 100. In this last case the two results are essentially undistinguishable in the figure. 

We can see that <5o(^o) displays first a strong decrease from zero down to large negative values after which it 
increases back to small negative values, and even becomes positive in the case of large mass ratios. For large 03 
approximate analytical results may be obtained for the location of the minimum and its value. They both scale as 

m-|- /mi . The essential qualitative feature of this result, namely the strong decrease of the phase shift, is directly 
linked to the large positive value of the scattering length. For example it is qualitatively fairly analogous to the one 
found for the I = phase shift of a square repulsive potential with strength Vq and range i?, which is given by: 



6o(k) = —kR + arctan 



k 



- tan(fc'-R) 



(42) 



where k is related to the energy E by E = h 2 k 2 /2m and E — Vq = h 2 k' 2 /2m. This is plotted in Fig(2] for 
{2mV a ) 1 ' 2 R/h = 10. 

We consider now the consequences of the above result for the s-wave scattering amplitude which is obtained from 
Eq.([9]). We are mostly interested in its modulus |/o(fco)| = | sin(5o(fco))| Ao- The result is plotted in Fig. [3] for 
m-f/rni — 100. The very interesting feature is the existence of several zeros occuring when SqQsq) = nir. The first 
few zeros are merely located at feg = m:/a 3 . Physically these zeros imply that there is no contribution from s-wave 
scattering to the total cross section. We see also that the scattering amplitude is strongly peaked for zero energy, a 




feature which is even more important for the corresponding contribution <7o(fco) = 47ra 2 /Q(fco) to the cross section. 
This feature is clearly linked to the strong decrease of the phase shift when ko starts from zero. As mentionned above 
this is directly linked to the large value of the scattering length, and consequently this is clearly not an artefact of 
our approximate treatment. This behaviour implies that a rough model for the s-wave component of the effective 
fermion-dimer interaction could be a constant up to the energy (hit /2as) 2 /2//y followed by no interaction at all for 
higher relative kinetic energy. This is quite different from the simple picture of a constant interaction, whatever the 
energy, which is valid for example when the fermions have equal masses. 

These results we have obtained appear quite naturally in the small m^/m^ limit. Their qualitative features are 
interesting enough to check that they are not just a curiosity of this simple, but not so physical, limit and that they 
rather stay valid qualitatively in the physical range for atomic mass ratios. We have checked that our result is in 
semiquantitative agreement with the exact numerical result from Eq.([5|. This is displayed first for m^/ra^ = 100 
which is somewhat beyond realistic values. We see from Fig. [3] that our approximation is in fair agreement with the 
exact one for reasonably small energies (the discrepancy at and near zero energy is just due to the difference between 
our approximate value for the scattering length and the exact one). There is however an increasing disagreement for 
higher values of the energy, which is not so surprising. Our approximate phase shift is varying with energy somewhat 
less rapidly than the exact one. 

If we go now to the more realistic value corresponding to a hypothetical mixture of 6 Li and 173 Yb displaying a 
Feshbach resonance, we obtain for this mass ratio m^/ra^ = 173/6 ~ 29 the results reported in Fig. [4] It is quite 
interesting to see that the results are not deeply modified compared to the low energy part of the preceding figure. 
Both the approximate and the exact result display a single zero. Just as in the preceding case, the exact zero is at 
a somewhat lower energy than what our approximate analytical expression gives. It is interesting that the relative 
location of this zero is in this order. Indeed our figure shows that the approximate zero is near the dissociation 
threshold of the dimer, and that it is no longer present for slightly lower values of the mass ratio m^/m^. On the 
other hand the exact zero is only halfway to the dissociation threshold. 

Since the existence of this zero is a landmark of a strongly varying s-wave scattering amplitude, it is of high interest 
that it exists for even lower, and accordingly more realistic, mass ratio m^/m^. We have found specifically that the 
exact zero is disappearing for m^/m^ ~ 15.27. If we take the lighter fermion to be 6 Li, this would imply for the 




FIG. 3: Modulus of the I = scattering amplitude |/o(&o)| as a function of ko, for mass ratio m^/raj. = 100. Dashed line: 
asymptotic result. Full line: exact numerical result. 



heavier one a mass quite near the one of Rubidium. Unfortunately the two stable isotopes of Rubidium are bosonic, 
but there are two fermionic ones, 86 Rb and 84 Rb, which have half-life of order of a month. Finally if we take the case 
of the 6 Li - 40 K mixture, which is displayed in Fig. [SJ we see that the scattering cross section (obtained by squaring 
the amplitude) near the dissociation threshold is only 5% of its zero energy value. 

Finally it is worth noting that the disappearance of this zero in the scattering amplitude when the mass ratio is 
lowered is strikingly similar to the same disappearance [BJ[H] in the function a 3 (k, 0), which comes in the calculation of 
the scattering length 0.3. In this last case the disappearance occurs at a somewhat lower value of the mass ratio, and 
0,3 (k, 0) has no simple physical meaning, which makes the interpretation difficult. However one could go continuously 
from this function to the scattering amplitude by varying the last parameter from to k. Hence the two functions 
are closely related, and the simple physical explanation for the zero in the scattering amplitude (it occurs when the 
phase shift reaches tt) gives a corresponding understanding for 03 (k, 0). 



III. P - WAVE SCATTERING 



We turn now to the p-wave contribution to the scattering amplitude. The analysis is done with the Skorniakov 
and Ter-Martirosian equations as detailed in section II A The only difference is that, instead of averaging Eq.Q to 
obtain the s-wave component, in order to obtain the p-wave component a3 P (fc, ko) — (1/Att) J aX^P^k, k )a 3 (k, k ), 
where -Pi(k, ko) = k- ko/Zcfco is the t = 1 Legendre polynomial, we have to project on this I = 1 Legendre polynomial. 

The resulting equation is quite similar to the s-wave one Eq.((5j). This is seen quite clearly by introducing the 
notations (making again use of reduced variables, but omitting systematically the bar over quantities in reduced 
units): 



Ct(k, q) = 1 - Rk 2 + k 2 + q 2 
/3(k,q) = R'kq 



(43) 




and defining the kernels: 



K s {k,q) 



1 



2/3(k,q) 



In 



ot(k,q) + fi(k,q) 



a(k,q) - f3(k,q) 



and 



K p (k,q) = 



1 a(k,q) 



In 



a{k,q)+ fi(h,q) 



a(k,q) - (3(k,q) 



(44) 



(45) 



Then the equations for the s-wave and p-wave components can be written: 



I? 



a3s, P (fc, k ) 



2 



dq 



K s , p (k, k ) 
q 2 a3s. v (l'k ) 



«, 



i8 



K s , P (k,q) 



(46) 



where a3 S (fc, ko) = 03 (k) used in section [TT| 

Since we are interested in this part in mass ratios smaller than in the preceding section [TTJ we provide first the results 
for the s-wave component for the specific mass ratios of interest. The results for |/o(fco)| in terms of ko are given in 
Figj6j for the mass ratios m^/m^ = 1, 2, 4, 8. We have also inserted the result for m^/m^ = 6.64 corresponding to the 
6 Li - 40 K mixture, already displayed in Fig. [5| Except for the fairly slow increase of the atom-dimer scattering length 
with the mass ratio found at k = 0, all the results are of order unity for k = 1 and keep decreasing for increasing k . 

We turn now to the p-wave component results. We first give in Fig[7|the results for mass ratios m^/raj, = 1, 2, 4 and 
6.64, stopping at the value corresponding to the 6 Li - 40 K mixture. For equal masses the p-wave component is always 
small, as it could be expected, and it can clearly be omitted when one deals with the atom-dimer scattering properties. 
Hence the atom-dimer vertex can, to a large extent, be taken as a constant proportional to the atom-dimer scattering 
length. No complication is thus expected to arise from this side when many-body properties will be investigated. 
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The situation is qualitatively similar for m^/m^ — 2 and even 4: the p-wave component is fairly small and can be 
neglected compared to the s-wave contribution. However, we see that this p-wave component rises steeply when the 
mass ratio goes from 4 to 6.64. This last ratio is at the border of the resonant domain. 

This domain is now displayed in Fig. [8j where the results for mass ratios m^/m^ = 6.64, 7., 7.5, 7.75, 7.875 and 8. 
are plotted. We see that the hump, present for m^/m^ = 6.64 around fco ~ 0.45, develops into a strong resonance 
at lower and lower energy when m^/m^ is increased. Although we have sampled mass ratios which are very close to 
each other, we see that the resonance grows very rapidly with increasing mass ratio. The physical origin is clearly 
the development of a virtual bound state at positive energy. The lifetime of this state is directly related to the width 
of the resonance. Since the resonance peak gets quickly narrower, the lifetime grows very rapidly with increasing 
mass. The link with the bound state which appears [3] at zero energy for m^/m^ — 8.17 is confirmed if we look at 
the position of the resonance peak as a function of the mass ratio. This is shown in Figj9] We see that the position 
A)q of the resonance peak extrapolates to zero for a mass ratio which is in close vicinity of 8.17. More precisely when 
the resonance peak reaches &q = 0, it will be infinitely sharp, corresponding to an infinite value for |/i(fco = 0)|. 
Hence in this case a^ p (k, 0) will be infinitely large, which implies that the homogeneous Eq.(46) (i.e. without the term 
K p (k, 0) in the right-hand side) has a solution for zero energy ko = 0. This is just stating the well-known result that 
bound states are solutions of the homogeneous integral equation corresponding to Eq.(46). It is easy to find the lowest 



ra^jra^ for which this homogeneous integral equation has a solution. We find that this occurs for m^/m^ = 8.172 in 
full agreement with Kartavtsev and Malykh [5]. 

Let us now come to the total scattering cross section a(ko), which is likely to be the easiest direct physical quantity 
to measure experimentally. Since we neglect angular momenta higher than I = 1, we have for this cross section: 



g(fco) 

47T 



= ^(2^ + l)|/,(fc )| 2 = l/o(fco)| 2 + 3|/ 1 (fc )| 5 



(47) 



The result is plotted in Figfl0| for ra^jra^ — 1,2,4 and 6.64. We see that, while up to m^/m^ = 4, the cross section 
is almost featureless, it display a strong resonance for m^/m^ — 6.64 which corresponds to the experimental value for 
6 Li - 40 K mixture. Naturally when the mass ratio is further increased, this resonance becomes even stronger. This is 
shown in Fig. 11 where we display also the case of mf/mj, = 8. The resulting resonance dwarfs the preceding one. 
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FIG. 6: Modulus of the 1 = scattering amplitude |/o(&o)| as a function of fco, for mass ratios m-f/mj. = 1, 2, 4, 6.64 and 8. The 
ratio m^/m_i = 6.64 corresponds to the 6 Li - 40 K mixture. All these quantities are in reduced units, equivalent to set a = 1. 



This strong resonance in the atom-dimer scattering cross section for 6 Li - 40 K mixtures should imply a number 
of experimentally observable consequences. The more natural way to evidence it would be to measure the collision 
properties between a 40 K cloud and a 40 K - 6 Li dimer cloud, in a way analogous to the recent experiments on 6 Li 
for spin states [30 . However this might not be so easy to perform. Measuring the equation of state of a mixture 
is on the other hand a fundamental question which should be easier to answer. The resonance should affect in an 
important way the low temperature equation of state of a mixture of 40 K cloud and a 40 K - 6 Li dimer, arising from 
an unbalanced mixture of 40 K and 6 Li atoms on the BEC side of the Feshbach resonance. This will not appear at 
the mean field level since the angular average will cancel the effect because this is a p-wave resonance. But going to 
next order in 40 K, analogously to second order perturbation theory, there will be no cancellation. Hence one should 
see a strong dependence on the K density. One expects also a strong density dependence to arise when the Fermi 
wavevector of the 40 K Fermi sea reaches the wavevector corresponding to the resonance, since for lower wavevectors 
the scattering is essentially negligible. One should see something similar to a threshold effect in density. We note 
that the resonance should also have an important effect even on a balanced mixture of 40 K and 6 Li, on the BEC 
side, when the temperature is raised. Indeed the dimers will be partially broken by thermal excitation which provides 
a natural source of free 40 K atoms, and effects analogous to the ones arising in the unbalanced mixture. Similarly, 
going to high temperature, we expect this resonance to have a marked effect on the virial coefficients since these are 
systematically related to the three-body problem [3T] . Naturally the resonance is also expected to affect strongly the 
transport properties. Both the effects on the equation of state and on the transport properties should appear in the 
frequency and the damping of the collective modes pQ. 

Another way to understand the effect of this resonance on the 6 Li - 40 K mixtures is to remark that the long lived 
virtual bound states responsible for it should behave in many respects in a way analogous to real bound states, as it 
has already been noted by Levinsen et al [23] . Hence on the BEC side the physical description should involve not only 
free fermions and dimers, but also the existence of trimers. Obviously all the physical properties should be affected 
by the presence of this additional fermion species. A related point is that the existence of these trimers should affect 
the dimer-dimer scattering properties. In other words we have only explored the simpler fcrmion-dimer scattering 
properties, but we expect that the dimer-dimer scattering amplitude will display related energy structure, and they 
should not be so complicated to explore with the methods of our preceding work [3]. 
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FIG. 7: Modulus of the 1 = 1 scattering amplitude |/i(fco)| as a function of fco, for mass ratios m,f/m± = 1,2,4 and 6.64. The 
ratio rn^/m_i = 6.64 corresponds to the 6 Li - 40 K mixture. 



Finally it is worth stressing that the use of optical lattices should provide a very convenient and powerful way to 
explore the resonance we have pointed out. Indeed the natural mass ratio of 6.64 between K and 6 Li atoms happens 
to be just at the border of the strongly resonating mass ratio domain. A rather weak optical lattice could be used 
to slightly tune the effective mass of 40 K or 6 Li (by making an appropriate choice of the light frequency, only one 
atomic species is affected) , producing strong modification of the resonance and hence of the physical properties of the 
mixture. In this way one could go into the strongly resonating domain, or even reach the threshold for the appearance 
of real bound states. Or one could go in the other direction and basically get rid of the resonance, which would 
allow to prove that it is responsible for specific physical properties of the mixture. More specifically the 3D potential 
V op t(r) = sEr \sm (Kqx) + sin 2 (/<o2/) + sin 2 (i^o- z )] with Er = i^ 2 /2rn provides, within second order perturbation 
theory, an effective m* given by m/m* — 1 — s 2 /32. Increasing the 40 K effective mass to reach the bound state ratio 
8, 17 requires s — 2.45 which corresponds to a fairly weak optical potential. This simple picture of effective mass 
modification works only if the size of the involved objects (dimer, trimer) is large compared to the optical wavelength 
A. Taking A ~ 500 nm this should be a valid approximation in the vicinity of the Feshbach resonance, where the 
scattering length becomes quite large. Even if corrections to this simple picture are necessary, the qualitative physical 
trends should remain valid. 



IV. CONCLUSION 



In this paper we have studied the s-wave and p-wave contributions to the fermion-dimer scattering amplitude as a 
function of the mass ratio. When masses are equal the physical situation is quite simple. The p-wave contribution is 
completely negligible, and the modulus of the s-wave contribution is essentially constant. Hence the energy dependence 
of the scattering properties is inessential, and the fermion-dimer scattering length a 3 is enough to fully characterize 
these properties. 

When the mass ratio is increased the situation becomes much more complex, both for s-wave and for p-wave. For 
the s-wave contribution we have obtained an analytical solution in the asymptotic limit of very large mass ratio. This 
allows to have explicitly the behaviour of the s-wave scattering amplitude. We find that it displays a large number 
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of zeros. These are directly linked to the very rapid increase of the phase shift with energy, a zero appearing each 
time the phase shift crosses a multiple of n. At low energy this rapid increase of the phase shift is directly linked to 
the known large value of the fcrmion-dimer scattering length for large mass ratio, although we find that at higher 
energy the phase shift saturates and goes back toward zero. In the range of mass ratio corresponding to possible 
experimental situations, the numerical solutions display a behaviour in qualitative agreement with the asymptotic 
limit, although the agreement is not quantitative which is not surprising. Interestingly we find that a zero is still 
present in the s-wave scattering amplitude down to an experimentally reasonable mass ratio of m^/m^ ~ 15. 

For the p-wave contribution the behaviour gets rapidly quite complex as soon as the mass ratio has a sizeable 
value. This is basically due to the appearance of a fermion-dimer bound state when the mass ratio reaches the value 
m-f/mi — 8.17, and we have limited our study to mass ratios below this threshold. Even below this threshold the 
existence of the bound state appears through the existence of virtual bound states at positive energy. These give rise 
to a resonance in the scattering amplitude. This resonance gets stronger and goes to lower energy when the mass 
ratio is increased toward the limiting value of m^/m^ = 8.17, at which the resonance corresponds to a divergence at 
zero energy. It is very interesting that the mass ratio m^/ra^ = 6.64, corresponding to the 40 K - 6 Li mixtures, much 
studied experimentally, is at the border where the p-wave resonance becomes important. Roughly, below this mass 
ratio we find a fairly small bump in the modulus of the p-wave scattering amplitude, which is fairly unimportant in 
the overall scattering properties and may be omitted in a first approach. On the other hand above this mass ratio, the 
p-wave resonance grows very rapidly and becomes the dominant feature of the scattering. We have stressed explicitly 
that the use of optical lattices could allow to vary experimentally the mass ratio throughout this very interesting 
range of mass ratio. One could in this way disentangle the role of this p-wave resonance in the physical properties 
of these 40 K - 6 Li mixtures, which are likely to be quite complex. In the same way one could check if these virtual 
bound states play an important role in the instabilities of these mixtures. 
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